procedure qr1k(input, output, nsum_in, offset)

string input
string output
int nsum_in
int offset

begin

string input_section
int size

nsum = nsum_in
imgets(input,"i_naxis2")
size = int(imgets.value)

input_section = input//"[*,"//offset//":"//size//"]"

print(input_section)

scopy (input_section,
output, w1=INDEF, w2=INDEF, apertures="1x3", bands="", beams="", apmodulus=0,
format="multispec", renumber=no, offset=0, clobber=no, merge=no, rebin=yes,
verbose=no)

end




